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Abstract 

A generalized version of standard map is quantized as a model of quantum chaos. It is shown 
that, in hyperbolic chaotic regime, second moment of quantum level velocity is ~ 1/h as predicted 
by the random matrix theory. 
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I. INTRODUCTION 



Chaotic systems are characterized by positive lyapunov exponents such that in phase 
space neighbourhood trajectories diverge exponentially. Among the chaotic systems there 
are special class of systems which are completely chaotic or hyperbolic |]| . Phase space of 
hyperbolic systems has only unstable orbits as there are expanding and contracting real 
directions with positive and negative lyapunov exponents respectively. Area preserving 
maps like cat map and baker map are examples of hyperbolic systems. Even the well known 
standard map of kicked rotor, a text book paradigm of Hamiltonian chaos Q, is not proven 
to be hyperbolic even for strong external kick strength. 

In the study of chaotic quantum systems, it is of fundamental interest to characterize 
highly chaotic (but not known to be hyperbolic) and hyperbolic chaotic regimes in quantum 
domain. This forms our motivation here, and to pursue further it would be more appropriate 
to quantize a single dynamical system which has parameters for highly chaotic and hyperbolic 
chaotic regimes. In one of our earlier works, a generalized version of standard map was 
introduced to study the dynamics of a kicked particle trapped inside an one dimensional 
infinite square well potential • The generalization, arising from length scales namely well 
width and field wave length, has parameters to fulfill the present requirement. With our 
knowledge, there is no other single system possessing parameters for the above mentioned 
classical regimes. 

In quantum domain, dynamics of levels in parameter space is known to have manifesta- 
tions of classical complexity. While quantum levels cross each other for regular case, they 
exhibit avoided crossings when underlying classical dynamics is chaotic. Level dynamics 
can be described by level velocity wherein system parameter plays the role of pseudo-time. 
In Ref. |i| the notion of curvature i.e., second derivative of levels with parameter, is in- 
troduced to quantify the avoided crossings. It is known that the curvature distribution of 
chaotic system follows an universal behaviour [Hf. 

In this connection, one another quantity of importance is the second moment of level 
velocity. A semiclassical analysis of kicked system shows that second moment is the sum 
of all classical time correlations of the kicking potential, such that lowest (zeroth) order 
correlation being the Random Matrix Theory (RMT) predicted second moment [6j]. It is 
also known for quantized standard map that, in highly chaotic regime there are systematic 
deviations between second moment and the corresponding RMT prediction. The deviations 
are thus related to the non- vanishing higher order classical time correlations 

In the present work, we introduce quantum version of a generalized standard map as a 
model of quantum chaos to study the level velocity statistics. In particular, we compute 
the second moment and compare with the RMT prediction in order to characterize different 
chaotic regimes. 

II. THE MODEL 

A. Classical system 

Considering a particle trapped in a one dimensional infinite square well potential V Q (q) 
of unit width (hard walls at q = ±1/2), which experiences a periodic kick from an external 
pulsed field. The Hamiltonian is 
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H = \ + %(,) + ^ (i) 

and underlying kick to kick dynamics of the particle is equivalent to discrete dynamics 
described by a dimensionless area-preserving mapping: 

Pj+i = Pj + 7^ sin(2vrrg j ) 

Qj+i = 9j + Pj+i ( 2 ) 

which is defined on 2-torus i.e., a unit square [—1/2, 1/2) x [—1/2, 1/2) with periodic bound- 
aries. Here r = 1/A is ratio of two length scales of the system namely, well width and field 
wavelength; k is effective strength of the kick. This is the Generalized Standard Map (GSM) 
which was introduced in our earlier studies on the above Hamiltonian 0]. 

GSM is continuous when r is integer and discontinuous otherwise. One can immediately 
recognize that widely studied standard map of kicked rotor is a special case (r = 1) of GSM. 
Since the standard map is a continuous map, for small (k < 1) dynamics is predominantly 
regular wherein many rotational invariant circles (also called as KAM tori) are interspersed 
in the phase space. They act as forbidden barriers for chaotic orbits to diffuse. Gradual 
destruction of these invariant structures with the increase of k, leads to onset of chaos; for 
k ^> 1, dynamics is highly chaotic. On the other hand, when r is non- integer no KAM tori 
exist in the phase space. In this case, depends of the parameter r the phase space is either 
mixed or fully chaotic even for small k values. 

The Jacobian J of GSM is such that 

I Trace J| = |2 + kr cos(2irrqj) \ . (3) 

Since \qj\ < 1/2, for r < 1/2 |Trace J| > 2. That is to say Jacobian has real eigenvalues. In 
other words, the system is completely chaotic or hyperbolic for r < 1/2. In this regime there 
are contracting and expanding real directions or alternatively stable and unstable manifolds 
throughout the phase space. Thus GSM is realized as a rare class of dynamical system as it 
has parameters for both highly chaotic and hyperbolic chaotic regimes. 



B. Quantum system 

GSM arises from the equation of motion of free particle in presence of a field V(q) which is 
applied as time periodic impulse. The field is defined as: V(q) = /ccos(27rrg)/(47r 2 r); V(q) = 
V(q + 1) and the Hamiltonian is 

H = ^+V(q)^6(j-t). (4) 

j 

By integrating Shrodinger equation over unit time we obtain corresponding quantum prop- 
agator as 

(J = e-^l^ e~ iV ^' n . (5) 
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Then the quantum dynamics can be described as \ip(t + 1)) = U\tp(t)), which is a quantum 
analogue of the classical map. 

On quantizing 2-torus phase space by introducing periodic boundary conditions both in 
q and p [Tf we have: q\n) = (n/N)\n) ; p\m) = {m/N)\m) where n,m = —N/2, —N/2 + 
1 . . . N/2 — 1. Here iV = (2irh)~ 1 is the dimensionality of the Hilbert space and the semi- 
classical limit is N — > oo. The position and momentum eigenstates obey the periodicity 
\n + N) — \n) ; \m + N) = \m) and the transformation function is 



(n\m) 



exp 



i2nmn 
N 



(6) 



Being a homogeneous linear system, Shrodinger equation in finite N dimensional space has 
solutions \4>j), j = 1,2, ...N, which are linearly independent. Since the Hamiltonian is 
time periodic (unit period), according to Floquet theory the solutions satisfy eigenvalue 
equation U\<f>j) = e~ % ^ \</>j) . Eigenstates \(f>j) are quasienergy states and eigenangles <fij 
are quasienergies. Then general solution at a given time is = Ylj c j\4>j) where Cj are 
constants. As a consequence of hermiticity of Hamiltonian, quasienergy states are orthogonal 
and they form a complete set in finite N dimensional space. 

Matrix form of the propagator in discrete position representation is 



U„ 



(n\U\ri) 



exp 




The Hamiltonian (J3J has reflection symmetry about the origin i.e., H(q,p) 
This symmetry is reflected in the quantum propagator matrix 



Un.n.' 



2nr 



N 



cos 



2nr 



(n f + a) 



through the relation [U, Ft] — where the hermitian operator Ft is defined as 

R\n) — | — n) for a = 
= |— n — 1) for a = 0.5 . 



(7) 

H{-q, -p). 
(8) 

(9) 



Note that a phase factor a is introduced in the matrix element to avoid exact quantum 
symmetry. Since Ft 2 = 1 we may label the eigenstates of U with eigenvalues ±1 of Ft i.e., 
the states are \<f>±). For a = 0.5, symmetry matrix of order N is 



R N = { n \Ft\n') = 5(n + n' + 1) (mod N) 



(10) 



which has ones along secondary diagonal and zeros elsewhere. Then the state components 
have a relation (— n — l\(b) = ±(n\d>) i.e., 



=> 




(11) 



Eigenstates can be numerically obtained by diagonalizing the matrix U nn > of order N. If 
is even integer, there are N/2 even parity states and N/2 odd parity states 

On exploiting -B-symm etry the diagonalization can be reduced to the matrix of order N/2 

by standard procedure (lOj. The reduced matrix is 
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e^/ 4 [-ikN \2ttv / 1 N\l) 
U nn i = — 1=- exp < cos \n H > 

x ^(n-n'f/N _j_ e in(n+n'+lf/N^ ^) 

where n, n' = 0, 1, ... iV/2 — 1. Now the separation of parity states is obvious. 

III. SPECTRAL STATISTICS 

One of the standard statistical measures for a chaotic quantum system is the nearest 
neighbour spacing distribution of quantum levels. For a regular system, levels are clustered 
such that the spacings follow Poisson distribution. On the other hand, in chaotic case 
the levels exhibit repulsion such that spacings exhibit RMT predicted Wigner distribution. 
We expect the classical complexity of GSM, arises due to the parameter r, will also have 
manifestation in the spectral spacings. It is evident from Fig. ^ that, in contrast to the 
predominantly regular case, for r = 0.5 (hyperbolic) the spacings follow Wigner distribution. 




FIG. 1: Nearest neighbour spacing distributions of 1000 quasienergies which correspond to even 
parity states for k = 0.3 with a = 0.5. Smooth curves are the Poisson and Wigner distribution. 
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A. Level velocities 



Having seen the effect of r in the spectral spacing, we further investigate on dynamics of 
the quasienergies in parameter space i.e., level velocities. To be specific, we study second 
moments of level velocities in different classical regimes. We take a = 0.35 so that R- 
symmetry is broken in the quantum system, and the factor a is dropped out in following 
expressions for the sake of convenience. Quasienergies (fij = (f)j(k,r) have scaled velocities: 



2ir 2 r 2 \ 


90, 


N ) 


dk 


27rV\ 




N ) 





= nr cos(2irrn/N) | (n|0 7 -) | 



(13) 



and 



Vj 



/2vrr 2 \ d(pj 



Nk 



\ Nk J dr 
{2nr 2 \ rl 



h 



{4>j\dV/dr\ 



-2nr J2( n / N ) sin(27trn/N) \ (n\<j>j) | 2 - J2 cos(27rrn/iV) | (n\ 



(14) 



Average velocities in semiclassical limit are 



(x) 

(y) 



— J2 X J= sin(Trr) 

3 

I sr^ t \ 2sin(7rr) 

— > Vj — cos 7rr 



(15) 



Then the second moment of x is given by 
1 



<* 2 > = ^E*J 



(vrr) 2 |]Tcos 2 (27rrn/iV) (|(n| 



+ Yl cos(2nrn/N)cos(2irrn'/N) (| (n]^) | 2 | (n'|0i> ' 2 



(16) 
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Assuming that spectral averaged eigenfunction components are independent of specific po- 
sition eigenvalues n, terms within the angle brackets can be taken out of the sum. Then 



(x 2 ) ~ (Trr) 2 [(|(n|0,.}| 4 )-(|(n|0 i }i 2 |(n / |^-}| 2 )]Ecos 2 (27rrn/iV) 



+ (|(n|0 i }i 2 |(n'|0 i }| 



cos(2irrn/N) 



(17) 



In chaotic regimes, standard RMT results 

(|(n|0,)| 4 ) = 3[N(N + 2)]- 1 ^3N~ 



(18) 



which correspond to Gaussian orthogonal ensemble are applicable here as well. It should 
be noted that application of RMT results essentially adopt the assumption made above. 
Replacing the sum by integration in semiclassical limit we arrive to 



( X )rmt 

Similarly the second moment of y is 



7rr 



N 



1 + 



sin(27rr) 
2nr 



+ <*> 2 • 



(19) 



(y 2 ) = (20) 

3 

= J2{i 27Tr ( n / N ) sin(2vrrn/iV)] 2 + cos 2 (2vrrn/A0 

n 

+ 4irr(n/N) sm(2Tcrn/N) cos(2vrrn/iV)} (\{n\4>j)\ l 



+ J2 {(27rr/N) 2 nn'sm(2irrn/N) sm(2irrn/N) + cos(2vrrn/iV) cos(27rrn' / TV) 

+ 4nr(n/N)sm(2nrn/N) cos(2nrn f /N)} ({(nl^flin'l^)] 2 ) (21) 
and the RMT prediction is 



if (nr) 2 sin(27rr) 



(y>BMx = ^i + - 3-+ 27rr 



5 / 

2 - OH 



\2 



3 



2 cos(27rr) + (y) 2 . (22) 



In chaotic regime, quantum states are such that the quantities in left hand side of Eqn. 
()18|) fluctuate about the respective RMT values. These fluctuations could lead to the failure 
of RMT in predicting the second moment. In Ref. a semiclassical analysis on the 
systematic deviation between the second moment and its RMT prediction has been made. 
Here we are interested to see the validity of the RMT prediction in different classical regimes. 
For that we define normalized deviations: 
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(x 2 ) - (A 

(x 2 ) 



A, 



(y 2 ) - (y 2 ) r 
(y 2 ) 



and taking average of the two positive quantities as 



A 



A x + A % 



(23) 



(24) 



The deviation is calculated for various parameters and plotted in Fig. |21 Let us first 
consider the data obtained for small values of k (0.3 and 0.1). As we see the deviation is 
nearly zero for r < 1/2. For further values of r, the deviation is found to be significantly 
large. This is expected as the underlying classical dynamics is predominantly regular or 
mixed and the RMT is not applicable. On the other hand, the data for large k (5 and 
25) exhibit different behaviour. Though the classical phase space is highly chaotic for these 
parameters, significant deviation observed for various values of r shows the failure of RMT 
prediction. However, it is noticable from the four sets of observation that for r < 1/2 the 
deviation is nearly zero. 




FIG. 2: Deviation between level velocity second moments and RMT predictions. In all the cases 
we have taken a = 0.35 and N = 200. 

In Fig. 01 the deviation is plotted for two cases by varying k. For r = 1 the deviation 
has strong fluctuations with k. The large deviation for k < 5 is due to the predominantly 
regular /mixed behaviour of the underlying classical system. We also observe significant 
deviation when k is close to integer multiples of 2tt. This may be attributed to the presence 
of "accelerator modes" [l^ . which are small regular regions embedded in the sea of chaotic 
phase space. On the other hand, for r = 0.5 the deviation remains negligible irrespective of k 
values. This shows that RMT predicts second moment of level velocity in hyperbolic chaotic 
regime. From the Eqn. (fT^j) . we see that (x 2 ) RMT ~ 1/N or ({d<j)/dk) 2 ) KMT ~ iV ~ 1/h. We 
obtain the same result for the velocity with respect to r as well. 
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FIG. 3: Deviation between level velocity second moments and RMT predictions. We have taken 
a = 0.35 and N = 200. 

B. Sawtooth map 

In this section we study another quantum model whose classical counterpart has hyper- 
bolic regime. Consider a free particle that is subjected to a time periodic impulsive potential 
V(q) = — Ag 2 /2; V(q + 1) = V(q). Kick to kick dynamics of such a particle is described by 
the sawtooth map [lj : 



(25) 



Pj+i = Pj + M 3 

defined on a unit torus. This map is stable for —4 < A < and unstable (or hyperbolic) 
otherwise. Quantized version of this map can be obtained, as usual, upon introducing 
periodic boundaries in q and p. Details of quantized sawtooth map are given elsewhere 
jl4j |. Note that here also the earlier symmetry arguments hold. From the corresponding 
quasienergies and quasienergy states, we define the scaled level velocity as 



Z 3 



NnJ d\ 



dV 



dX 



Y.(n/Nf\(n\^)[ 



(26) 



with average (z) = 1/12. Then the second moment is given by 



(^ 2 ) = T?J2 z j=H HN) 4 (|H<^>| 4 ) + £ (n/iV) 2 (nViV) 2 (Kn|0,)| 2 |(n'|0,)| 2 ) . (27) 



Repeating the earlier procedures, we find the RMT prediction as 
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FIG. 4: Deviation A' for the quantum sawtooth map with a = 0.35 and N = 200. 



Defining normalized deviation as 



(A 



A' 



40N 



+ <*> 5 



(z 2 ) - (A 



(28) 



(29) 



in Fig. 0]we have plotted the deviation for different A values. In the stable region —4 < A < 0, 
where the RMT result is not applicable, the deviation is large. Evidently, RMT predicts 
level velocity second moment in hyperbolic regime. Eqn. (J28J) shows that {z 2 ) 

rmt ~ 1/iV or 

((d<j)/d\) 2 ) RMT ~ N ~ 1/h. Here also we see that second moment is inversely proportional 
to the Planck constant for hyperbolic chaos. 

IV. CONCLUSION 



In this paper we have introduced quantum version of a generalized standard map. The 
classical system corresponds to the dynamics of a particle, trapped in an ID infinite square 
well potential, in presence of time-periodic kicks. For different classical regimes, the second 
moment of quantum level velocity is computed and compared with the RMT prediction. 
We have shown that while the prediction fails in highly chaotic regime, the second moment 
is well predicted by the RMT as ~ 1/H in hyperbolic chaotic regime. By considering 
another example viz., quantum sawtooth map, the RMT prediction of second moment in 
hyperbolic chaotic regime has been reinforced. We hope the presented result would shed 
new lights to explore more on quantum hyperbolic chaos. 
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